Testing models of inflation with CMB non-gaussianity 
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Two different predictions for the primordial curvature fluctuation bispectrum are compared 
through their effects on the Cosmic Microwave Background temperature fluctuations. The first 
has a local form described by a single parameter /jvl- The second is based on a prediction from 
the warm inflationary scenario, with a different dependence on wavenumber and a parameter fwi- 
New expressions are obtained for the angular bispectra of the temperature fluctuations and for the 
estimators used to determine /atl and fwi- The standard deviation of the estimators in an ideal 
experiment is roughly 5 times larger for fwi than for /jvl- Using 3 year WMAP data gives limits 
—375 < fwi < 36.8, but there is a possibility of detecting a signal for fwi from the Planck satellite. 



I. INTRODUCTION 



Inflationary models of the early universe have proved to be in spectacularly good agreement with observations 
of the Cosmic Microwave Background. The spectrum and the gaussianity of the temperature fluctuations are both 
consistent with the simplest type of inflationary models. The quality of the observational data will improve in coming 
years, and a future aim will be to learn more about the nature of the inflationary era. 
f~i , We shall focus here on the possible use of non-gaussianity to distinguish between different inflationary scenarios. 
C^' The dissadvantage with this approach is that the amount of non-gaussianity produced in the simplest inflationary 
Q models is very small ^j, \2i, lS] , and it might not be possible to separate it from non-primordial sources of non-gaussianity 
$_i > such as background sources Q or the second order Sachs- Wolfe effect [1, Attention has therefore concentrated 
^ [ on inflationary models which produce a relatively large amount of non-gaussianity and, fortunately, there are many 
. examples. 

An early example was the curvaton mechanism, where the density fluctuations originate from the quantum fluctu- 
ations of a light scalar field whose energy density is negligible during inflation 0, S| . The non-gaussianity depends on 
^ ' the energy density of the curvaton when it decays, and it may be larger than secondary sources of non-gaussianity 
r in the CMB in special cases [1]. Another possibility is that non-gaussianity can be produced if there exists a non- 

\l homogeneous reheating mechanism after inflation, and the non-gaussianity may be comparable to secondary sources 
. The possibility which we focus on here is motivated by a recent result on warm inflation (To| , in which it is possi- 
' ^ ■ ble to have models where the primordial non-gaussianity is significantly larger than the non-gaussianity produced by 
C ' secondary sources. 

The non-gaussianity can be quantified by taking moments of the density fluctuations, especially the bispectrum 
B^{ki, fc2, /ca). The bispectrum is defined by the stochastic average of a triple product of perturbation amplitudes 



>: (C(ki)C(k2)C(k3)) = (2^fSc(fci,fc2,A;3)<53(ki+k2 + k3). (1) 

"-H . 

' Note that we will use the Bardeen variable C, (the three-space curvature on a constant density hypersurface [ll|) to 
JH \ describe the primordial density fluctuations. The Bardeen variable has the advantage of being constant for small 
■ - - ■ amplitude density perturbations on large length scales. 

The shape of the bispectrum depends on the particular inflationary scenario. In this paper we will concentrate on 
two different models for the bispectrum. The first will be the local model, 

^(fci, fe, A;3) = 2/c/ Pdki)Pdk2). (2) 

eye 

where 'eye' indicates a cyclic permutation of fci, ^2 and and the primordial power spectrum P^(fc) is defined by 

(C(ki)C(k2)) - (27r)3F^(fci)<5=^(ki + k2). (3) 

This type of bispectrum describes the primordial fluctuations for inflationary models with density fluctuations gen- 
erated by a curvaton field jl2| . The usual parameter which is used in this case is not fci but the non-linearity 
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parameter /at^ of the Newtonian potential fluctuations. The best observational limits on the non-linearity at present 
comes from from the WMAP three-year data [13i] , 

- 54 < Jnl f /c/ < 114 at 95% C.L. (4) 

The Planck satellite observations have a predicted sensitivity limit of around IJnl] ~ 5 [1]. 

The second model is motivated by the non-gaussianity arising from the warm inflationary scenario [13] ■ In warm 
inflation there is a significant amount of radiation present during inflation and the density fluctuations arise from 
thermal, rather than quantum fluctuations (l4l. [TH- The bispectrum takes the form 

B^' iki,k2,k3) = fwiJ2iK^ + K^)Pdki)Pdk2)l^i -^2. (5) 

eye 

A similar bispectrum appeared in pf^ as part of a longer expression describing second order perturbative effects 
following the usual inflationary scenario, but there fwi ~ 3, and here we consider arbitrary values of fwi- We shall 
consider the prospects for testing this model using CMB data from WMAP and PLANCK experiments. For other 
forms of bispectrum, see refs. \lLd8,M- 

II. THE ANGULAR BISPECTRUM 

We follow ref. (l6j and express the primordial bispectrum as a Legendre polynomial expansion, 

Bc(fci,fe,A;3) = 5]^/i(fci,fe)Pc(fci)^c(fc2)Pi(ki -ka). (6) 

eye I 

Suppose that the /;(fci, fc2) are homogeneous functions of momenta, then without any loss of generality, we can write 

B^ki, k^M) = /"'^"'(^i' ^2, fcs), (7) 

where 

eye 

The CMB temperature fluctuations have an angle-averaged bispectrum Bi^i^i^, defined by 

^ I mi m2 <«'i™i«'^™^«'3™3> (9) 



where the matrix is the Wigner 3j symbol and temperatme fluctuations have been expanded in spherical harmonics, 

^Y.^imYUi^). (10) 

The coefficients in the spherical harmonic expansion are related to the curvature fluctuations by 

^C(k)5/(k)y4(k) (11) 

(Note that many references on the microwave background use the transfer function from the Newtonian potential 
fluctuation. We use the curvature fluctuation, which is currently used by CMB FAST software [20].) The angle- 
averaged bispectrum for the curvature bispectrum ([7]) will decompose into a series of terms, 

Bhhl3 = ^ fniBll\^i^ (12) 

n.l 
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where, from eqs. ([J), ([9]) and (fTTI) . 



TO2 W3 



f_,-Vl+i2+/3 J_ Ml 
771177127713 ^ 

/ n ('^^fc*5i,(fci)ii:77i,(k7)) B"'(fci,fc2,fc3)<5(ki +k2 +k3) (13) 

i=i 

This can be simphfied by expanding the delta and Legendre functions in terms of spherical harmonics and eliminating 
the momenta. After using Wigner 3 — j symbol identities (see ref [2l|), the result is that |16j 



/•OC 

perm;';' " 



(14) 



where 'perm' denotes a permutation of h, I2 and I3. The dependence on the transfer function is contained in two new 
functions 

bl], - - / dkk'^+'Pc{k)gi{k)j,{kr) (15) 
5f ^ - - / dkk^gi{k)ji{kr), (16) 

where our notation is based on ref (Most recent papers have used Greek letters for the 6i's, e.g. ai for bf^, but 
this notation does not extend very well due a shortage of usable Greek letters.) 

The coefficients c{l'il'2l; hhh) were first obtained by Komatsu They are related to Wigner 6j symbols, 



(2;i + l)(2Z2 + l)(2Z3 + l), 

(17) 



Z'l /3\ Mi ^1 M ^2 '2 ^ \ f/2 ^1 ^3 

J \ o oj \ o oy\/i / 

The 6j symbols vanish unless the l'^ lie in the range hi I, and therefore this expansion is most useful for small values 
of I. In such situations, eq. (|14p combined with a good Wigner 6j evaluation code can become an efficient way to 
evaluate the angle-averaged bispectrum. 

There is a useful normalisation condition on the coefficients, 



'1*2 



A. Local models 



The local bispectrum corresponds to the model with I = n = Q, which we will denote by the label 'CP. In this case 
I'i = k, i = 1)2, and the coefficients c{l[l'2l; hhh) reduce to 

The angle-averaged bispectrum is simply 

POO 

Bflis = E ^lil^O-Jihk) I bZ\r)bZ\r)br^\rydr, (20) 

perm ^ 

where we define 

2 

bf\r)^- dkeP^{k)gi{k)ji{kr). (21) 

Jo 
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The functions 6f^^(r) and 6p^(r) are evaluated using CMBFAST or related software [20|- Some care has to be 
exercised in the evaluation of (r) . It is important that the time evolution of the density perturbations in CMBFAST 
starts early in order to catch the large k behaviour of the transfer function correctly. The function, bf^^ (r) has a large 
spike not far from r = , where 



ait) 



dt 



(22) 



is the speed of light times the conformal time difference between decoupling td and the present t^. This is shown in 
figure [T] The resolution in r has to be very high to capture this accurately, but this incurs a large time penalty when 
evaluating the angle averaged bispectrum for many different values of I. We get around this problem by evaluating 
bf^^{r) with a high resolution (3200 values of r in the range ztrd) and then averaging to produce a function bf^{r) 
with a lower resolution. The integral over r can be obtained accurately using bf^{r) because the other terms in the 



integrand vary comparatively slowly with r. 
integral are described in ref [17|] 



More sophisticated methods of reducing the number of points in the 
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FIG. 1: The function fef ^(r) is plotted as a function of r in Mpc for I — 2, I = 100 and / = 250. The cosmological model is the 
WMAP ACDM plus n = 1 model, which has nA = 0.788, = 0.212 and h = 0.778. With these parameters, r* = 1.41 Gpc 
and Td = 323 Mpc. 



B. Non-local models 



The second model which we consider in detail is the model with I — n~ 1, which we will denote by the label 'WI'. 
The coefficients c{l'il'2^] hhh) can be found using explicit formulae for the 6j symbols, given in ref [2l|, for example. 
They can be expressed conveniently in terms of Edmond's angle 0, defined by 



cos 6* = 2 



/l(^l + l) + /2(^2 + l)-^3(/3 + l) 
(2Zi + l)(2/2 + l) 



(23) 



This angle is an interior angle of the simplex with sides corresponding to the entries in the 6j symbol. The non- 
vanishing coefficients are given by 



cih±l,l2±i,l;h,l2,h) 
cihhO-Jihh) 



1 

— I cos( 
4 



1 ± 



1 



± 



1 



1 



2h + l 2^2 + 1 (2/i + l)(2;2 + l) 



(24) 



5 

= -T COS0 - 1 T i^r-T ± l^r-T + ^o, , 1^.o, . 1^ ■ (25) 



cihhO-Jihk) 4V 2/i + l 2^2 + 1 (2^i + l)(2;2 + l) 

The terms depending on r can be simplified by using Bessel function identities, 

bW^ = l:b?'{r) - hf^'ir) (26) 
bU - ^6p^(r)+6r'(r) (27) 
bul, = Ur{r)-br'{r) (28) 

= ^6r(r)+5r'(0, (29) 

where the prime denotes a derivative with respect to r and bf'^{r) is given in eq. (|2ip . We have introduced a new 
function bf^^ (r), defined by 

2 

bT'{r)^- dkP^{k)gi{k)n{kr). (30) 
This is plotted in figure together with the angular power spectrum C/, defined by, 

Ci^- I dk^ P^{k)gi{kf. (31) 




FIG. 2: The Ci, bf'{r,) and are plotted as a function of I. The C; and hf'{r,) have been scaled by l{l + l)/(27rA2), 

where = 2.5 x 10"^, and the {r*) by an additional (/ - 1)(/ + 2)/r^ The cosmological model is the WMAP ACDM plus 
n = 1 model, which has Q.a = 0.788, fim = 0.212 and h — 0.778. With these parameters, r, = 1.41 Gpc and rd = 323 Mpc. 



The angle-averaged bispectrum Bj^i^^ is given by substituting the above functions into eq. (|14p , 



perm 

dr 



[(li + (^2 + feP^(r)6,^^(r)cos0 + 6,f '(r)C^'(r)r2) 6/^^(0, (32) 

where c{lil20;lil2l3) is defined in eq. bf^{r) is defined in eq and &p^(»') in eq. (PTjl . The radial hmctions 

and their derivatives can be evaluated using CMBFAST. 
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III. MODEL TESTING 
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FIG. 3: The theoretical standard deviation A/ for an ideal experiment is plotted as a function of the maximum value of I for 
the two models of non-gaussianity discussed in the text. The cosmological model is the WMAP ACDM plus n = 1 model, 
which has Qa = 0.788, = 0.212 and h = 0.778. 



We shall use Maximum likelyhood methods methods to consider the sensitivity of 'ideal' experiments for detecting 
the non-local model of the bispectrum which we introduced in the previous section, and we shall consider how well 
this model can be distinguished from the local model. For our purposes, an ideal experiment is one in which there is 
no noise in the detector and no secondary sources of non-gaussianity. 

For a single model, the log-likelyhood function ~X^{B°^ili^, f) is a function of the nongaussianity parameter / for 
a set of observations B?^,^, . 



For an ideal experiment, we take the simple anzatz 



E 



( oobs 



(33) 



where Ci is the angular power spectrum and Ci-^Ci^Ci^ is the theoretical variance of the bispectrum when the non- 
gaussianity is small. The sumations extend to a maximum / value Imax which represents the limiting resolution of the 
experiment. 

This form of likelyhood function is associated with an estimator 



6F ^ 



oobs D 



Ci^Ci^Ci, 



where the normalisation factor F is the Fisher information. 



F 



(34) 



(35) 



For an ideal experiment, the mean of / is the value of the parameter / and the variance of / is F^^ . 

The standard deviation A/ = has been plotted in figure[31 using the values of Bi^i^i^ discussed in the previous 

section. The WMAP ACDM plus n — 1 cosmological model has been used to define the cosmological parameters. The 
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standard deviation for the local model is in good agreement with previous calculations, even when they use different 
cosmological parameters [11] . 

The calculation of the standard deviation for the WI model with Imax = 512 takes about 30 minutes on a 2GHz 
Intel Core Duo processor. The time is taken up by the high resolution used to evaluate the radial integrals when we 
get to values of I ~ 400. The high resolution is needed mostly for the 'squeezed triangles', where one value of I in 
the summand is far smaller than the other two. This fact shows up clearly in figure IH where the plots show that 
the dependence on resolution dissappears when low multipoles have been omitted from the sum. Our resolution is 
adequate for the maximum value I = 300 in the next section. For the Planck experiment we would need to push the 
limit up to ? = 2000 or so. In this higher range, an increase in the standard deviation arising from dropping low 
multipoles may well be acceptable as a trade-off for having to perform the calculations with many values of r. For 
example, the resolution used here is already sufficient to obtain the variance of the estimator for the WI model from 
modes in the range 9 < Z < 1800, and it comes out to be Afwi = 15.9. This is an upper bound on the variance of 
the estimator which uses all of the modes 1 < I < 1800. 
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FIG. 4: The theoretical standard deviation A/ for the WI model is plotted as a function of the maximum value of I with 
different resolutions in the radial integrals and different lower limits on the values of I. The cosmological model is the WMAP 
ACDM plus n = 1 model, which has SIa = 0.788, ^2™ = 0.212 and h = 0.778. 



The observations can also be used to distinguish between two different models 



information is now a matrix. 



1,2,, 



6 ^ 

I1I2I3 



Ci^Ci^Ci^ 



The Fisher 



(36) 



The separability of two different models is associated with a low value for the correlation 

(This quantity was called the 'cosine between two distributions' in ref. '23*1). If two bispectra have good separability, 
then the failure to detect any non-gaussianity for one bispectrum does not rule out the possibility of detection of 
non-gaussianity in the other bispectrum. The square of the correlation for the CI and WI models has been plotted 
in figureini At / — 300, the square of the correlation is about 10%, which is quite low but not low enough to separate 
the two models definitively using this statistic. 
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FIG. 5: The square of the correlation r for an ideal experiment is plotted as a function of the maximum value of I for the two 
models of non-gaussianity discussed in the text. The cosmological model is the WMAP ACDM plus n = 1 model, which has 
Ha = 0.788, f^™ = 0.212 and h = 0.778. 



IV. NON-GAUSSIANITY TESTS ON CMB DATA 

The estimator for the non-gaussianity of model i in an ideal experiment was given in the previous section, 

-1 Dobs ryi 

?idcal ^ \p ^lihl-ihhh 

'lf2'3 

where the observed angle-averaged bispectrum is related to the spherical harmonic coefficients by 



DObs 



E 



h h h 

TOl 7712 TT^3 



^/imi ^l2va2 ^Ism^ • 



(38) 



(39) 



We can make use use the fact that the theoretical bispectrum eq. (pOI or eq. ((5^ consist of a small number of terms 
which factorise to devise an efficient method for evaluating the estimators based on ref. [23 |. 
The Wigner 3j symbol is first replaced by the integral 



h h I. 

TOl TO2 TT, 

where c{lil20;lil2l3) is defined in eq. to give 

f ideal 



^2 ) = cihhO-Jihh)-' I d^nYi,^,ih)Yi,^,{n)Yi,rnAi^), 



Now we define a radial sequence of filtered sky-maps cr(r, n) by 

Ci 



l.rn 



6^-f(r) 

awi{r.,-a) = E "St — a;m^m(n) 



l,m 



Ci 



(40) 
(41) 

(42) 
(43) 
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For the local model, Bf^J^i^ is given by eq. After interchanging the sum and the integrals in cq. ([^T|l . 

/cr' =^ I dPnr^dr Gci<JciOf,L (45) 



This can be evaluated easily with the HEALPix [26] software package for analysing CMB data. 

The second model has angle averaged bispectrum given by eq. (p2|) . The 1(1 + 1) terms in Edmond's angle Q ((23)) 
can be replaced by the angular Laplacian D^, defined by 

D'Yi^{xv)^~l{l + \)Yi^{xv). (46) 

The estimator becomes 

fwf = IT^— / c'^ndr {(p^aci)(ywiONL + o-c/(£'Vw/)crjVL - OdOwA^'^ <yNL) - 2aci(Jwi'^NLr^} ■ (47) 

We can regard this as a volume integral and use the divergence theorem, 

fwT' = I d^'x (Vac/) • {^cjwi)cJNL, (48) 

where V is the three dimensional gradient operator. This expression can be evaluated at high speed using the 
subroutines for handling angular derivatives of CMB data in the HEALPix software package. 



A. Real experiments 

Galactic and extragalactic sources contaminate the cosmological microwave signal and have to be masked out, 
leaving us with a set of spherical coefficients aim of the masked sky. The statistical properties of these coefficients are 
affected by the masking procedure. An important change appears in the theoretical covariance matrix 

Ci 

imil2m2 — (fliimiaigmg >, (49) 

which is no longer diagonal. Detector noise is also an important consideration. The detector noise is not uniform over 
the sky and the noise in the microwave intensity has a non-trivial covariance matrix. 

Creminelli et al. [13] have introduced an improved estimator which should be close to optimal for the CMB 
observations. Their estimator 

h = + (50) 

where the linear piece 

~ ~'2N Cl Cl Cl ^'l™l'2m2'^'3'"3- (''I) 

lll-llz 12 3 

The normalisation factor N = Fu/skyj where /sky is the fraction of the sky which remains after the mask is applied. 

The linear part of the estimator can be expressed in terms of filtered maps by following the same steps as in the 
beginning of this section. We shall concentrate here on the effects of the mask. To a reasonable approximation, it can 
be shown that 

blM.Yl,mAi^)Cl,ra,l2rn2Yi,raAi^)^M{n)Y,^-hjCl (52) 

where hi is any slowly varying function of I and the mask function M(n) = in the region being masked out and 1 
otherwise. The details are given in the appendix. Using this approximation, the linear parts of the estimators become 



1 

n'ci 



fl^i ^ -— I d^nr^drM{2SciNLCTci + SciciaNL) (53) 



fwi ^ j^Jd^n r'dr M (^Sciwi^nl + s'^^I^jCJnl + ^/iVL^W/ + S^winl^'ci) (54) 
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where the radial functions SAsir), A, B = CI, WI or NL, are defined by 

2l + lht{r)hf{r) 
4^ Ci 

2l + lbfirybf{r) 



(55) 
(56) 
(57) 
(58) 

The linear parts should be combined with the other terms pS]) and ([48]). 

Crenelli et al. [3 have argued that since the contributions from the linear terms depend on angular averages of the 
filtered maps, they can be reduced in size by arranging that the average value of the data outside the mask is zero. 
However, we have found that the contributions from the linear terms are important and they need to be included. 



Sab {t) 
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SAsir) 
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S^Alir) 


I 




-E 



4^ Ci 

2l + lbfirybf{rY 
,4. a ■ 



B. WMAP data 



We shall consider the non-gaussianity estimators for the 3-year data which has been made available by the WMAP 
collaboration. The angular resolution of this data should allow values of I up to 300-400, but we have chosen to take 
the maximum value of I for which detector noise can be neglected, and this is around / — 300. We have used 800 
values of r in the range 13750Mpc < r < 14350Mpc and 200 values of r for 8750Mpc < r < 13750Mpc. (Values of r 
outside this range appear to make very little change to the estimators). 

The values of the estimators and the variance of the ideal experiment are shown in table IH The rms values of 
the estimators for a set of randomly generated gaussian maps have been used to check that the linear term in the 
estimator removes the effects of the mask without biasing the estimator. Leaving out the linear term gives a spurious 
signal. 

The results for fci can be expressed as a constraint on the more commonly used parameter /tvl, where fci = 
0.6 Jnl, 

-26</jvL<109 at 95% C.L. (59) 

which we should stress ignores detector noise. This is consistent with the results published by the WMAP collaboration 
given the fact that the detector noise increases the standard deviation of the estimator by 10-20% [1^. The results 
for the WI model are 



- 375 < /w/ < 36.8 at 95% C.L. (60) 

Although there is no evidence for a signal, it may be worth pointing out that fwi < at 94% C.L., and this would 
be in agreement with the predictions of the warm inflationary model [lo| . 



Map Mask 


/sky 


I 


fci 


Afci 


fwi 


Afwi 


W kpO 


0.765 


300 


24.9 


20.5 


-169 


105 


V kpO 


0.765 


300 


33.7 


20.5 


-174 


105 


R kpO 


0.765 


300 


-0.712 


13.4 


3.09 


89.2 



TABLE I: Values of the estimator (with the linear correction) and variance for a selection of 3- year WMAP datasets and 
masks. The frequency bands are denoted by W and V and the maps are 'foreground reduced'. The row labelled R gives the 
average and rms average of the estimator over 20 randomly generated gaussian maps. The deviation A/ in rows W and V is 
the deviation for an ideal experiment scaled by the unmasked fraction of the sky fsky 
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V. CONCLUSION 

We have compared two different bispectra which represent the non-gaussianity in the CMB produced by two different 
types of infiation. The first bispectrum (CI) is of the local type, which can result from a curvaton field and the second 
(WI) was modelled on a result deduced from the theory of warm inflation. In the second example we found a new, 
simple form for the estimator used to quantify the amount of non-gaussianity present in the CMB data. 

The WI bispectrum demands accurate handling of the numerical integrations and also careful consideration of the 
effects introduced by masking out part of the sky. We have introduced a binning technique to improve the accuracy 
of the integrations and used an improved estimator to deal with the masked data. 

The variance of the estimator for the WI nongaussianity in an ideal experiment using modes up to ^ = 500 is plotted 
in figure m Beyond / ~ 500, we would have to improve the resolution in the integrations, remove some of the lower I 
modes or adopt an alternative computational strategy [13] ■ For example, using modes in the range 9 < ^ < 1800 sets 
an upper bound on the deviation of Afwi < 15.9. 

When applied to the WMAP 3-year data we reproduce similar results to previous work for the CI bispectrum. So 
far there is no evidence for the WI bispectrum, and the limits are 

- 375 < fwi < 36.8 at 95% C.L. (61) 

According to the calculation in ref. J2| : the strong regime of warm inflation leads to the WI bispectrum with 

/H./«-181n(l + ^) (62) 

where r is a large parameter. It is possible to predict the probability of a 2cr signal from the Planck satellite for this 
model. Using Afwi < 15.9, a 2a signal would be obtained with better than 95% probability if r > 65. 

APPENDIX A: SOME PROPERTIES OF MASKED DATA 

Consider the situation where we expect a gaussian signal (T(n) but due to unwanted background sources we mask 
out an area of the sky Ai . The masked spherical harmonic coefficients are 

ai,n^ j d^nM{n)a{n)Yil^{n) (Al) 

where M(n) = for n G and Af(n) = 1 otherwise. The covariance matrix of the masked coefficients is defined by 

This can be expressed in terms of the angular power spectrum Ci of the full-sky coefficients, 

CHm.i.m^ = E / rf'Aid2n2M(ni)M(n2)r,:„^(ni)r,;„^(n2)r/™(ni)r/™(n2). (A3) 

Im 

The main objective of this appendix is to approximate the sum 

Pi^)= KbhYhniA^)Chmihni:2YhniA^)^ (A4) 

/i/2inim2 

where bi is a slowly varying function of I. Using eq. (IA3I) . 

F{h)^ y"d2n,d2^2M(ni)M(n2)i^(n,rii)if(n,n2)^C/r/™(ni)y,™(ri2) (A5) 

Im 

where the kernel K{n. n') is given by 

if (n, n') = J2 biYi„An)Yi,^in'). (A6) 

Im 

The next step is to clarify the sense in which the 5/ are slowly varying. Introduce a new function b{t) by 

5/ = / e-'('+i)*5(i)dL (A7) 
Jo 
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We say that the bi are slowly varying if b{t) = for t > fi, where /i is a small parameter. For primordial fluctuations, 
the scale of variation of the foj's is set by the acoustic oscillations, and /x is a small parameter determined by the 
conformal time at the last scattering surface. (In reallity b{t) will not vanish for t > /i, but it must be relatively small 
in view of the smoothness of the plots of the fe;'s). 

The kernel K{n, n') is related to the fundamental solution to the heat equation on a unit sphere K{h, h' , t) by 



K{n,n') = / dtK{n,h' ,t)b{t). 
Jo 



For small time, we can use the asymptotic expansion 

i^(n,n',t)^^e-'^("'")^/4*(l + 0(t2)), 
where d(n, n') is the great circle distance between n and n', to get an upper bound, 



(A8) 



(A9) 



(AlO) 



where Kq — J^ii'^^ + 1)^' ^'^d 0{rf') corrections are dropped. 

If the mask is empty, then the orthogonality property of the spherical harmonics can be used to show that the 
function F takes a constant value S, where 



(All) 



Now use eq. (jA5|) to calculate the difference 

\F{h)- SM{n)\< J d^hid^h2 (M(ni)M(n2)-M(n))i^(n,ni)if(n,n2)^Gyi™(ni)l^,„(n2; 

Im 

We have a bound (|A10p on the kernel and we can bound the final sum by a constant C, 



(A12) 



Im 



< 



I 



Consider the case where n £ A4, which forces the integrand into the rest of the sky A4, 

1 



|F(n) - S'M(n)| < CKI [J rf2fj/g--d(n,n')V4r, j 



(A13) 



(A14) 



There is an assymptotic expansion for the integral for small r\ obtained by setting n' = Ac + e, where fie is the closest 
point on the edge of the mask to n. The result is that 



|i^(n) - 5^/(6)1 < CKt,^x{2^Ye 



(A15) 



A similar argument applies when n e A^. The approximation F k, SM (h) holds as long as the point on the sky does 
not approach too closely the edge of the mask where d{n, he) — 0. 
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